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Abstract 

The phenomenon of Anderson localization is studied for a class of one- 
particle Schrodinger operators with random Zeeman interactions. These 
operators arise as follows: Static spins are placed randomly on the sites of 
a simple cubic lattice according to a site percolation process with density 
x and coupled to one another ferromagnetically. Scattering of an electron 
in a conduction band at these spins is described by a random Zeeman 
interaction term that originates from indirect exchange. It is shown rig- 
orously that, for positive values of x below the percolation threshold, the 
spectrum of the one-electron Schrodinger operator near the band edges is 
dense pure-point, and the corresponding eigenfunctions are exponentially 
localized. 

Localization near the band edges persists in a weak external magnetic 
field, H, but disappears gradually, as H is increased. Our results lead us 
to predict the phenomenon of colossal (negative) magnetoresistance and 
the existence of a Mott transition, as H and/or x are increased. 

Our analysis is motivated directly by experimental results concerning 
the magnetic alloy EurCai-zBe. 



1 Introduction 



Theoretical understanding of the effects of disorder on electron transport 
in metals and semiconductors started with Anderson's famous 1958 paper 
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[3]. In this work, Anderson argued that in the presence of strong disor- 
der caused by impurities and/or defects and neglecting electron-electron 
interactions, electrons populating a weakly filled conduction band of a 
metal get trapped in exponentially sharply localized one-particle orbitals. 
This is the phenomenon known as "Anderson localization". A conse- 
quence of localization is that the conductivity of such a material very 
nearly vanishes at low temperatures. If disorder is described by an on-site 
random potential with bounded distribution and short-range correlations 
in a one-electron tight-binding Hamiltonian then Anderson's arguments 
can be made precise, mathematically, for one-dimensional systems with 
arbitrarily weak disorder [12, 15], and for higher-dimensional systems, 
provided that disorder is strong enough or the energy of the one-electron 
orbital lies in the band tails [10, 9]. It is generally expected — but not rig- 
orously proven — that, in two-dimensional systems of this kind, all states 
are localized, no matter how weak the disorder is. In contrast, in three 
or more dimensions, localized states with energies in the band tails are 
expected to coexist with extended states (generalized eigenstates of the 
model Hamiltonian) corresponding to energies in the continuous spec- 
trum near the center of the band, provided the disorder is sufficiently 
weak. It is expected that wave packets made from superpositions of such 
extended states exhibit diffusive propagation corresponding to a non-zero 
conductivity [21, 1]. One is led to predict that, at very low temperatures, 
a three-dimensional disordered semiconductor exhibits a transition from 
an insulating state (all electrons in the conduction band occupy local- 
ized states) to a conducting state (some fraction of the electrons populate 
extended states), as the density of electrons in the conduction band is 
enhanced or the strength of disorder is lowered. This transition from an 
insulator to a metal is called a "Mott transition". 

The purpose of this paper is to describe and analyze a tight-binding 
model for a Mott transition where the disorder is caused by indirect ex- 
change interactions between the electrons in a conduction band and a 
dilute array of localized atoms with non-vanishing spin. This model has 
been introduced in order to interpret theoretically electronic properties of 
Europium-based hexaborides (EuxCai-^Bg), [23]. A Mott transition is 
found experimentally as the concentration, x, of the magnetic Europium 
atom is varied. This feature is captured correctly in our model. Fur- 
thermore, the model also suggests an explanation of the phenomenon of 
"colossal negative magneto-resistance" encountered in experimental ex- 
plorations of EuicCai-ajBe; see [23]. 

Fairly simple arguments presented below lead us to introduce a model 
given in terms of a one-electron tight-binding Hamiltonian with a random 
Zeeman interaction term acting on electron spin. This term describes 
indirect exchange interactions between an electron in the conduction band 
and electrons in the half-filled 4/ shell of a Eu-atom located nearby. It 
takes the form of a ferromagnetic coupling of the spin of the electron in 
the conduction band to the static total spin of electrons in the 4/ shell of 
a Eu-atom. Because the latter is quite large, S = 7/2, it can be described, 
in good approximation, by a classical unit vector, in [16]. However, if a 
unit cell of the simple cubic lattice of a Eu^Cai-^Be alloy contains a Ca- 
atom then rh = 0, because a Ca-atom has spin 0. At low temperatures, 
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the direction of m is approximately constant throughout a connected Eu- 
cluster, because of indirect ferromagnetic exchange interactions between 
the spins of different Eu-atoms in the cluster. However, the direction of 
fn varies randomly from one Eu-cluster to another, as long as the external 
magnetic field vanishes. Thus, an electron in the conduction band of a 
Eu^Cai-zBe alloy in zero magnetic field propagates in a disordered quasi- 
static background of essentially classical spins located in those unit cells 
that contain a Eu-atom. These spins are ferromagnetically coupled to the 
spin operator of the electron. 

One of our main contentions in this paper is that, as long as there 
is no ferromagnetic long-range order (unit cells containing a Eu-atom do 
not percolate), but the concentration of Eu-atoms is not too small, in zero 
magnetic field, this type of magnetic disorder causes Anderson local- 
ization in the tails of the conduction band. 

If the concentration, x, of Eu-atoms is brought above the percolation 
threshold then there is an infinite connected cluster of positive density 
of unit cells containing a Eu-atom, and the alloy is observed to order 
ferromagnetically at low enough temperatures [24]. Most Eu-spins are 
then aligned in a fixed direction. The same happens if a sufficiently strong 
external magnetic field is applied. Finally, if x is very small most unit cells 
exhibit a vanishing spin, that is, the vector rh vanishes in most unit cells. 
In all these three situations, the disorder felt by electrons in the conduction 
band is weak, so that the localization threshold (or "mobility edge" ) moves 
towards the band edges. We thus expect to observe a derealization- or 
Mott transition to a conducting state, as x increases across x c , or if the 
external magnetic field is increased. 

It is not understood, at present, how to prove the existence of such 
a transition and analyze its characteristics, although, heuristically, it is 
fairly well understood. But the existence of Anderson localization in the 
band tails, for x below the percolation threshold, but non-zero, and for 
a sufficiently weak external magnetic field, can be proven rigorously. A 
sketch of our proof is the main message of our paper; technical details can 
be found in [6]. We have profited from previous results in [4, 19]. 

While the main mathematical results presented in this paper may not 
be particulary suprising, they concern examples of Anderson localization 
that have not previously been studied mathematically. 

Our discussion is summarized Figures 1 and 2. 

Our paper is organized as follows. In Section 2, we briefly recapitulate 
some experimental findings on Eu^Cai-zBe and sketch physical mecha- 
nisms that may be at the root of the observed phenomena. They lead 
us to propose, in Section 3, some idealized models that are expected to 
yield an adequate theoretical interpretation of the observed phenomena. 
In Section 4, we formulate our main mathematical results on our models 
and sketch how they may be used to interpret various features found in 
the experiments. In Section 5, brief outlines of proofs of our main results 
are presented; (for details, the reader is referred to [6]). 
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Figure 1: x denotes the concentration of Eu-atoms, H the value of a homogeneous ex- 
ternal magnetic field. The concentration of conduction electrons is assumed 
to be approximately constant. The shaded area corresponds to an insulating 
state; a Mott transition to a semi-metal is expected to be observed at its 
boundary. Rigorous results are known for a subset of the parameter values 
inside the shaded area. 




Figure 2: a denotes the conductivity of the alloy. The figure provides a qualitative 
plot of a. 



2 Summary of experimental results con- 
cerning EiLrCai ^Bg, and physical mecha- 
nisms 

We begin by recalling some essential properties of EuBe. This binary 
compound crystallizes in a simple cubic lattice. At the center of each unit 
cell of the crystal there is a divalent Eu atom, at every corner of a unit 
cell there is an octahedron of B-ions; see Figure 3 below. 

The 4/-shell of a Eu-atom is half filled, which, according to Hund's 
rule, implies that the total spin is s = 7/2. Electron transport is domi- 
nated by defect-state conduction, with a low concentration, n c , of around 
1CP 3 charge carriers per unit cell [25]. At low temperatures, EuB6 or- 
ders ferromagnetically at a Curie temperature Tc m 12 K, accompanied 
by a significant reduction of the resistivity, p, in the ordered phase [8]. 
The isostructural compound CaB6 is obtained by replacing Eu by iso- 
electronic but non-magnetic Ca, which leads to a further reduction of n c 
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Figure 3: Schematic unit cell of EuE>6; each cube corner is the centre of a Boron 
octahedron 



by an order of magnitude [23]. In the series Eu^Cai-^Be, Tc decreases 
monotonically with decreasing x, down to x ~ 0.3. At lower values of 
x, no onset of long-range magnetic order is observed. Instead spin-glass 
type features dominate the magnetic response at low temperatures [24]. 
For the simple cubic Eu sublattice, x c — 0.31 is the site percolation limit 
[24]. In the concentration range 0.2 < x < 0.3, significant localization 
and colossal magnetoresistance effects, such as shown in Figure 4, have 
been observed. For x = 0.27, the enhancement of the low-temperature 
resistivity by six orders of magnitude below 10 K may be quenched by 
rather moderate magnetic fields of the order of 1 T. Detailed investiga- 
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Figure 4: Main panel: p(T) of Eu^Cai-zBe for various values of x. The thin solid line 
for x = 0.27 is to guide the eye. Inset: Magnetoresistance of Euo.27Cao.73B6 
at low temperatures. 
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tions using selected-area electron diffraction patterns and high-resolution 
transmission electron-microscopy (HRTEM) have shown that also for large 
concentrations of Ca for Eu, the structural quality, that is, the perfect 
atomic arrangement in a simple cubic lattice is preserved and the disor- 
der is simply in the spins on the sites of the Eu clusters. Energy-filtered 
TEM reveals a phase separation into microscopically small Ca- and Eu- 
rich regions, respectively. This implies that the material is magnetically 
and electronically inhomogeneous [23]. 

Next, we sketch some ideas on a possible mechanism that may explain 
the long-range ferromagnetic order observed in EuB6 at temperatures be- 
low Tc- See also [18] for a similar discussion. The large size of the unit 
cells of EuB6 (as compared to the size of a Eu atom) and numerical simu- 
lations [14] suggest that ferromagnetic order is established through indi- 
rect exchange mediated by electrons in a somewhat less than half-filled 
valence band, with strong on-site Coulomb repulsion preventing double 
occupancy; see Figure 5. For a non-vanishing density of holes in the va- 
lence band [25] , the spins of the electrons in the valence band are expected 
to order ferromagnetically at very low temperatures. For the groundstate, 
this is a prediction of the Thouless-Nagaoka theorem [20, 17, 2]; (see also 
[11] for an analysis of ferromagnetism in the Hubbard model). Because 
of overlap of the orbitals of electrons in the valence band with those in 
the 4/-shells of Eu-atoms, the spin of a valence electron in a unit cell has 
a tendency of being "anti-parallel" to the total spin of the Eu-atom in 
the same unit cell, provided the temperature is low. Appealing to Hund's 
rule, this is seen to be a consequence of Pauli's exclusion principle and of 
the half-filling of the 4/-shell. Hopping processes of valence electrons into 
either an empty orbital of the 4/-shell of a Eu-atom or to an empty orbital 
of the valence band thus give rise to ferromagnetic order among the spins 
of the Eu-atoms and those of the valence electrons, the latter being "anti- 
parallel" to the spins of the Eu-atoms. Because the orbitals of conduction 



electrons overlap with those of valence electrons, there are exchange in- 
teractions between conduction- and valence electrons that, because of the 
Pauli principle, favor anti-ferromagnetic order between conduction- and 
valence electrons. Thus, the spins of conduction electrons have a tendency 
of being aligned with the spins of the Eu-atoms. We will describe this ten- 
dency by a Heisenberg term that couples the spin of a conduction electron 
in a unit cell ferromagnetically to the spin of the Eu-atom in the same unit 
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Figure 5: Spin ordering in EuB 6 . 



G 



cell. Since the spin of the Eu-atom is rather large (S = 7/2), we propose 
to describe it as a static classical spin, m. It would be of considerable 
interest to improve the theoretical understanding of ferromagnetism in a 
one-band Hubbard model coupled to a lattice of large localized spins. 

Our main concern, in this paper, is to provide some qualitative un- 
derstanding of electronic properties of Eu^Cai-^Be. Experimentally, a 
Mott transition from a metallic state to an insulator is observed, as the 
concentration, x, of Eu is lowered. Furthermore, for x < 0.3, when a ho- 
mogeneous external magnetic field is turned on, "colossal negative mag- 
netoresistance" is observed, [23]. These are the phenomena we wish to 
focus on in this paper. 

In our somewhat idealized theoretical description of Eu^Cai-^Be, we 
place Eu- and Ca-atoms at the centers of the unit cells of the simple 
cubic lattice Z 3 according to a site percolation process, with probability 
x to find a Eu-atom at a given site. The mechanism for ferromagnetic 
order through indirect exchange described above suggests that, within 
connected clusters of unit cells filled with Eu-atoms, the spins of the Eu- 
atoms are ferromagnetically ordered. Since different Eu-clusters are sep- 
arated by regions filled with non-magnetic Ca-atoms, one expects that 
the directions in which the spins of Eu-atoms are aligned vary randomly 
from one Eu-cluster to the next, as long as the external magnetic field 
vanishes (or is very small) . If there is no infinite cluster of Eu-atoms this 
introduces disorder, and, because the conduction electrons are scattered 
at the spins of the Eu-atoms, it enhances a tendency towards Anderson 
localization of conduction electrons. 

The threshold for the emergence of an infinite connected cluster in a 
site percolation process on 1? is x c ~ 0.31. For x above x c , one expects 
that there exists an infinite connected cluster of Eu-atoms. At low tem- 
perature, the spins of the Eu-atoms in the infinite cluster are all aligned, 
so that spin-disorder is weak. But if x is small there is an infinite cluster 
of non-magnetic Ca-atoms, while Eu-clusters are tiny, on average, and 
sparse. Hence, spin-disorder is again weak. However, for x in some range 
below x c , and in zero external magnetic field, there is considerable disor- 
der in the way spins in different Eu-clusters are aligned. This enhances 
scattering of conduction electrons at different Eu-clusters, and one expects 
that the mobility edge, E,, separating low-lying localized orbitals from 
extended states near the center of the conduction band is shifted away 
from the band edge towards the center of the band. If the Fermi energy 
in the conduction band is approximately constant as x varies one is led to 
predict that Mott transitions may be observed at some x* > x c and some 

<§C x c ; see Figure 5. 

3 Idealized model for Mott transitions 
driven by spin disorder 

In this section, we propose a model expected to exhibit some of the phe- 
nomena described in the last section, namely the Mott transition and the 
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Figure 6: Mobility edge, E*(x), as a function of x. 



colossal magnetoresistance observed in Eu^Cai-^Be alloys. Our model is 
idealized to an extent that some of its properties, in particular Anderson 
localization, can be established rigorously. 

Because, experimentally, the conduction band of Eu^Cai-zBe is only 
weakly populated, n c < O(10~ 3 ), we neglect interactions among con- 
duction electrons and describe the propagation of a conduction electron 
with the help of a one-particle model. It is convenient to make use of a 
tight-binding approximation. The Hilbert space of pure state vectors of a 
conduction electron is then given by 

H = l 2 (Z 3 )®C 2 . (1) 

Although valence electrons mediate an indirect exchange interaction be- 
tween conduction electrons and the electrons in the half-filled 4/ shells of 
Eu atoms, they do not appear explicitly in our model. Instead, the inter- 
actions of conduction electrons with the local Eu spins are described by a 
Heisenberg term coupling the spin of a conduction electron to the spin of 
a Eu atom localized in the same unit cell. Since the latter is quite large 
(s=7/2), we describe it by a classical unit vector, m. The Heisenberg 
term then takes the form of a Zeeman term, — Jrh-3, where 3 is the vector 
of Pauli matrices associated with a conduction electron, and J > is a 
constant. If a unit cell j £ 1? is filled with a Eu atom then \rhj\ = 1; if it 
is filled with a Ca-atom then \ fhj\ = 0. Eu- and Ca-atoms are assumed to 
be distributed over the unit cells of Z by a site percolation process, 
with probability x to place a Eu atom at any given site. The configu- 
ration, {rn,j} JgZ 3, of classical spins is treated as quenched (in particular 
time-independent). Because of the observed tendency of Eu-spins in a 
connected Eu-cluster to order ferromagnetically, the distribution of the 
configurations {fhj}j &I fl of Eu-spins in every connected Eu-cluster, C, is 
chosen to be given by a Gibbs measure 

dPc(m) := Zc l exv{n ^ m, • rrij + f3H ^ mf} 5{\rhj | 2 - l)d 3 mj . 
i,jec j jec 

\i-j\=l 

(2) 

Here, k = k(T) is a temperature-dependent, positive constant (k(T) is 
decreasing with T), /3 is proportional to the inverse temperature, H is the 
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(3) 

strength of a uniform external magnetic field in the z-direction, m - is the 
^-component of fhj , and Zc (the cluster partition function) is chosen such 
that dPc is a probability measure. 

The distribution of the Eu-clusters, C, is given by an independent site 
percolation process with density x. 

One might envisage to combine the distribution of the Eu-clusters C 
and of the configurations {mj}j 6Z 3 (with fhj = if j is occupied by a 
Ca atom) into a single probability distribution that would then describe a 
tendency towards Eu-Ca phase segregation. To simplify matters, we will 
not consider this possibility in the present paper. 

The one-particle tight-binding Hamiltonian is chosen to be given by 

h(uj) :— T + Wj (cj) — Jrhj (cj) • a , (3) 

where T e = —A is a short-range hopping term (A is the discrete Lapla- 
cian), uj denotes the randomness of the interaction terms, {fhj(uj)} is 
distributed according to (2), and w{ui) is a Bernoulli random potential 
with distribution 

, . [ w if fhj ^ (that is, j occupied by Eu) 
Wj((jj) = < ^ (4) 

I — w if fhj = (that is, j occupied by Ca) 

The potential w is incorporated in (3) because the potential energy of a 
conduction electron at a site j may depend on whether j is occupied by 
a Eu atom or a Ca atom. 

The quantity of main interest to us, in this paper, is the electrical 
conductivity 

(5) 

where D is the diffusion constant of conduction electrons. At temperature 
T = and for a given Fermi energy Z?f, D is given by 



D= f F dEp(E)D(E) . 

J —OO 



(6) 



where p(E) is the density of states. In linear response theory, D(E) is 
given by the Kubo formula 

p(E)D(E) = \™^^2\j\ 2 E\{0\(h(u J )-E-is)- 1 \j)\ 2 , (7) 

where E denotes an expectation with respect to the distributions given 
in (2) and (4), with Eu-clusters constructed according to an independent 
site percolation process of density x. 

In order to keep our exposition simple, we consider some limiting 
regimes of the model introduced in (2)-(4). 

(A) k -J- 0, H small. 

This regime is appropriate to describe electronic properties of 
Eu^Cai^Bg in the absence of magnetic order (e.g. for x very small, or 
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well above the Curie temperature of the magnetic transition). Mathe- 
matically, this is the easiest regime. Using methods developed in [4], it 
is not very difficult to establish Anderson localization, for a fixed value 
of x > 0, provided the energy lies sufficiently close to the band edges 
(depending on x and the value of the constant J in (3)). Coexistence of 
localized states corresponding to energies in the band tails and extended 
states corresponding to energies near the center of the conduction band 
is expected for J small enough. However, the nature of the spectrum of 
the random Schrodinger Hamiltonian h(uj) defined in (3) near the center 
of the energy band is very poorly understood, at present. 

(B) k — ^ oo, H small. 

In this regime, the spins of Eu atoms in every connected Eu-cluster are 
completely aligned, but their direction can vary arbitrarily from one such 
cluster to another one. Anderson localization can be proven for energies 
sufficiently close to the band edges, for sufficiently small values of the Eu 
concentration, x, so that, in particular, Eu does not percolate. For x suffi- 
ciently close to I an infinite, ferromagnetically ordered Eu-cluster of den- 
sity fairly close to x exists, and we expect to find two mobility edges close 
to the band edges. (A mobility edge separates energies corresponding to 
localized states from energies corresponding to extended states.) Mathe- 
matically, the existence of mobility edges remains, however, an open issue. 

(C) H — > co. 

In this limit, all the spins rhj are aligned in the positive ^-direction. 
The conduction band then splits into two independent subbands for elec- 
trons with spin in the negative z-direction and those with spin in the 
positive z-direction, respectively. Within each subband, the Hamiltonian 
h(uj) is then equivalent to a "Bernoulli Hamiltonian" 



Adapting methods developed in [19], we show that, at all energies in small 
intervals attached to the band edges, except possibly in subsets of those 
intervals of very small Lebesgue measure, the quantity p(E)D(E) intro- 
duced in (7) vanishes, as long as x =fc 0, 1. However, it is not known 
whether corresponding eigenstates are exponentially localized. It should 
be pointed out that the localization effect in this regime seems to be very 
weak, as can be seen in the inset of Figure 4. 

Assuming that mobility edges, E*, exist, separating energies E with 
p(E)D(E) — from energies E' closer to the center of the band, where 
p(E')D(E') > 0, we expect (on the basis of our mathematical analysis 
of regimes (A) and (B)) that, as a function of H, E, = E*(H) moves 
ever closer to a band edge, as H increases (that is, as magnetic disorder 
decreases). Thus, for < x < 0.3 and for a small, but positive density of 



h(w) 



T + v(lo), 



(8) 



where 
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conduction electrons, it can be expected that, at zero temperature, our 
model describes a Mott transition from an insulating state at small values 
of the magnetic field H to a conducting state at large values of H . If 
correct this conjecture would explain the colossal (negative) magnetore- 
sistance observed at x = 0.27 and very low temperatures, recall Figure 
4. 

In the remaining sections of this paper, we switch gears, from physical 
reasoning to mathematics. We state some mathematical results proven 
for the model introduced in (2)-(4), in regimes (A),(B), and (C), and we 
present outlines of proofs. 

4 Main mathematical results 

For the model of the electronic structure of Eu^Cai-^Be introduced in 
Section 3, we are able to prove the following main results. 

In regime (A): 

Theorem 1. The spectrum ofh{ui) close to the band edges is almost surely 
pure point, and the eigenf unctions decay exponentially fast. A fortiori, 
p(E)D(E) vanishes for E in the pure point spectrum. 

In regime (B): 

Theorem 2. For x < x c , the spectrum of h(u>) close to the band edges is 
almost surely pure point, and the eigenfunctions decay exponentially fast. 
A fortiori, p(E)D(E) vanishes for E in the pure point spectrum. 

If the Fermi energy lies in the pure point spectrum, the alloy is an 
insulator. The region where we can prove pure point spectrum is shown 
to shrink with increasing x (decreasing disorder), suggesting the conjec- 
ture that conduction increases as the absolutely continuous part of the 
spectrum approaches the Fermi energy. 

In regime (C), we consider, as discussed, the Bernoulli Hamiltonian 
(for each subband) 

h(u>) = -A + Vj(u), (10) 
where we choose the energy scale such that 




v with prob x 

—v with prob 1 — x . 



Theorem 3. In the band tail from [—v, —v 2 ], p(E)D(E) vanishes for all 
energies outside a set of the order of exp(— exp(l / ^/v)) . 

This is a perturbative result, and is meaningful only for small values 
of v. 
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5 Outline of proofs 



Although rigorous proofs of the theorems stated in the previous section 
are somewhat beyond the scope of this article (see [4, 6] for details), we 
try to convey the main ideas underlying these proofs in the following. 

The resolvent, or Green function, G(E) = (h— E)^ 1 , is closely related 
to the unitary time evolution, e _1 , via the Laplace transform 

G(z)ip = i / e itz e- [th i>dt, 
Jo 

for Im 2 > 0. It is therefore not surprising that detailed knowledge of the 
resolvent is a means for investigating transport properties of the system 
described by the Hamiltonian h. We consider the random Hamiltonian 
h(oj) given in (3), but for clarity of presentation we will treat the two 
random terms separately, that is, we first set w = 0. 
The Kubo formula 

p(E)D(E) = Km^ J2 W 2 E|<0|(/*H-£-ierV>| 2 (11) 

XG1.3 

shows that in order to establish absence of diffusion, we have to control 
E|G(£ + iO;0,a;)| 2 , 

with 

G{z;0,x) := <0|(fc(w) -*)">>■ 

The way we do this is to prove that, for a suitable choice of the energy 
E, depending on disorder, \G(E + iO; 0, x)\ decays exponentially for large 
\x\, with probability approaching 1, as \x\ — > oo. This will enable us to 
deduce localization (that is, pure point spectrum of h, almost surely) in 
certain energy regimes. 

It is technically convenient to study a regularized version of the full 
Green function: We restrict the Hamiltonian to finite-size cubes A C Z 3 
(with appropriate boundary conditions), with the intention of sending 
A /■ I? in the end. That is, we consider the matrix Hamiltonians 

h A , = -A A , + l Al (j)(wj -Jrhj-ff), 

where Aa, is the finite-difference Laplacian restricted to A\ C 1? , a cube 
of side length I, and 1a { is the characteristic function of A t . We then 
define 

G Al (z) := (h Al - z)~ L . 

The two main techniques used in our proofs are perturbation theory 
around "good" cubes, using the second resolvent equation, and induction 
on the scale of the cubes — a so-called "multiscale analysis". A "good" 
cube is defined as follows. 
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Definition. A cube Ai is called "good" at energy E if the resolvent 
Ga 1 (E) is well-behaved, that is, if there exists a constant c > such that 

||G^(£)||<e il/2 (12) 
\G Al (E;x,y)\ <e- c l^l for \x - y\ > ± . (13) 

If a cube is not "good", it is called "bad". Condition (12) — a so-called 
Wegner-type estimate — ensures that E does not lie too close to an eigen- 
value of the matrix \ia u and condition (13) represents the off-diagonal 
decay typical for the resolvent of a Schrodinger operator for E in a spec- 
tral gap. 



5.1 Multiscale analysis 

The main idea of multiscale analysis is simple. The various scales are 
the sidelengths / of the cubes A\. We establish existence of "good" cubes 
with high probability at an initial scale, and use this information to show 
existence of "good" cubes with even higher probability, inductively, at 
larger scales. The probability that a cube Ai is "good" should behave 
like 1 — pi, where pi oc I , and k is some integer. The key point is 
that, at finite scales, we have to prove existence of "good" cubes only 
with high probability and not with probability one: Denote by Q the 
whole probability space, that is, the set of all configurations oj of the 
random magnetic moments {jri,-} JgZ 3 . From the outset, we can discard 
configurations u that are difficult to handle, as long as they have small 
probability. 

Induction step. The induction step from scale l n to l n +i < l n proceeds 
as follows: Consider a cube Ai n+1 of size l„+i. Since "bad" cubes of size l n 
are, by the induction hypothesis, improbable, the probability that there 
are many, say N, "bad" i n -cubes in a Z n+ i-cube is even smaller, of the 
order of p,^. The integer N is chosen such that pf n < pi n+1 . Therefore, 
only configurations ui where the l n +i-cube contains not too many "bad" 
l n -cubes need to be considered. 

Next, we describe the basic perturbation step. Consider a set X C Z d 
and Y C X. The second resolvent identity yields 

Gx = Gy © Gx\y + Gy © G x\y^Gx , 

where V is the boundary operator linking Y to X\Y . For example, for 
x € Y and y G X\Y , we have that 

G x (E;x,y)= GY(E;x,z)Gx(E;z',y). 

(z,z>)edY 

To get an estimate on \Ga, i (E;x,y)\ we proceed in the following way. 
In a first step, we excise the "bad" Z n -cubes (of which there are not too 
many) from Ai n+1 . Iteration of the resolvent identity along a sequence of 
Zn-cubes the first of which is centered at x, the second on the boundary of 
the first, and so on until the last contains y, and the induction hypothesis 
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(13) at scale l n establish the estimates (12) and (13) at scale l n +i for 
Ai n+1 , with the "bad" cubes excised. The difficulty is now that, to couple 
a "bad" cube A\ n back to Ai n+1 , we need a slightly larger cube covering 
Ai n that satisfies a Wegner-type estimate (12). But, by the induction 
hypothesis, the Wegner estimate holds for cubes of this size only with 
probability 1 — pi n , whereas we should like to establish it with the much 
larger probability 1 —pi n+1 - 

For random potentials with a bounded probability density, it has been 
known for a long time [22] how to establish a Wegner-type estimate (12) 

,1/2 

with probability 1 — e simultaneously on all scales I. After the incep- 
tion of the rigorous study of Anderson localization [10], much effort has 
been devoted to the study of Wegner estimates, and there have been many 
advances, particularly in the continuum (see [13] and references therein). 
On the lattice however, for more general probability distributions, or, as in 
our case, matrix-valued random potentials, Wegner estimates are difficult 
to establish. Recent mathematical work [4] , triggered by our study of the 
hexaboride alloys, shows how to establish a Wegner estimate inductively. 
Inductive Wegner estimate. We have already argued that we can 
restrict our attention to configurations uj where there are not too many 
"bad" cubes of size l n in Ai n+1 — call these "bad" cubes Bk- The key idea 
is to modify each configuration uj by changing in ui — {rhj}ieA l i only 
the values of frij for j inside the such that, for this new configuration 
uj' , each "bad" cube has a neighbourhood satisfying a Wegner estimate. 
In a second step, one shows with the help of complex analysis that these 
configurations uj' have actually very large probability. 

Technically, this is done as follows (see also Figure 7): Pick one of 
the "bad" cubes, call it B. By the induction hypothesis, we know that a 
neighbourhood of a "bad" cube B satisfies a Wegner estimate with prob- 
ability 1 — pi n (the fact that such a neighbourhood is actually slightly 
larger than l n is a technical nicety and of minor importance). Now, cover 
the "bad" cube B with a slightly larger cube -B(i). Then the probability 
that the configuration in the ring B^\B makes the cube B^ "bad", no 
matter what the configuration inside B, is smaller than pi n . Thus, the 
probability that there are many, say N, equicentered cubes B^ of increas- 
ing diameter such that, for all i, the configuration in the ring Br i - ) \B( i _i- ) 
makes the cube B^ "bad" , no matter what the configuration in the inte- 
rior, is very small, namely of the order of pf n . We can therefore, 
from the outset, restrict considerations to configurations u> where each bad 
cube B can be replaced by a larger cube B( io _ 1 ) with the property that 
the configuration uj can be modified inside the cube B^ io _ l - ) alone such 
that the cube B(j _i) has a "good" neighborhood B(i )- The modified 
configuration shall be denoted by uj' . 

We can now use again the resolvent identity to establish the bound (12) 
on \\Ga, +i (E,uj')\\. We have thus found that ||Ga, l (S)|| is bounded 
at scale l n +i for one fixed configuration uj' . One may think that this is far 
too little, since a single point ui' has zero probability in fl. However, in a 
last step, one shows, using a matrix-valued Cartan-type lemma, that the 
probability of configurations satisfying the estimate (12) at scale l n +i is 
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Figure 7: Generic cube of size l n+ i with not too many "bad" cubes of size I 
very large. 

To understand the last step we need to recall a result from complex 
analysis, known as Cartan's lemma. The precise mathematical statement 
can be found in the appendix. For our purposes, a rough understanding 
of the lemma will suffice: The lemma says that an analytic function that 
is bounded away from zero at one point of its domain of definition is at 
most points not too close to zero. We need a higher-dimensional matrix- 
valued analog of this lemma due to Bourgain [4] , which we relegate to the 
appendix. This generalized lemma says that if ||Gyi n+1 (E; u)\\ is not too 
large for one configuration u', it is not too large for most configurations; 
actually only exponentially few (in l n ) configurations have to be excluded. 
It is important to point out that in order to apply this lemma the distri- 
bution of magnetic moments on the unit sphere needs to have a bounded 
density with respect to the uniform measure on the sphere. 

We can now apply the lemma, for we have explicitly constructed a 
configuration u' with the desired properties. 

Using (12), the off-diagonal decay (13) is easily established as before 
by expanding the resolvent along a sequence of nested cubes. 

The induction step is thus completed, because we have verified that 
(12) and (13) hold with probability at least 1 — Pi n+1 at scale l n +i- 

Induction basis. To complete the proof we need to verify conditions 
(12) and (13) at some initial scale lo- The proof is based on the follow- 
ing intuition. If a configuration cj is to be such that Ga 1() (E, lo) has an 
eigenvalue close to the upper edge of the spectrum Eo = id + |m|, then, 
first, most of the sites in Ai g have to be occupied by a magnetic moment, 
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and, second, most of nearest-neighbor moments have to be closely aligned, 
to maximize the hopping term. But since the moments are independent, 
a large-deviations estimate shows that such a configuration occurs with 
very small probability. (Here we use the fact that the probability of large 
clusters is exponentially small in the cluster size, for x < x c .) 
Absence of diffusion. From the exponential decay of Ga,(E), with 
probability ~ 1 — const l~ k , we get absence of diffusion, using that 

lim Ga{E; x, y) = G(E; x, y) , 

and the Kubo formula (7) for the diffusion constant D. 
Pure point spectrum. Call I the set of energies where we can prove 
almost sure exponential decay of G(E; x, y) (that is, / lies in the band 
tails). To prove pure point spectrum in /, we need additional ideas, see 
[9, 5]. The first is that, with respect to the spectral measure of h, almost 
every energy is a generalized eigenvalue, that is, a polynomially bounded 
solution to the equation 

(h-E)tP = 0. 

Next, we note that we can express such a solution by the Green function, 
1>(x)= Yl G Rn (E-x,z)iP(z') Vx€R n , 

where R n is most conveniently chosen a ring, say A; 5n \/l; n . The expo- 
nential decay of the Green function in R n implies thus exponential decay 
of tl>(x) for x large enough, which, in turn, means that ip is a true eigen- 
function, and, a fortiori, that E is an eigenvalue. 



5.2 Bernoulli random variables 

If the random potential is of Bernoulli-type, the above inductive method 
for proving a Wegner estimate breaks down, since the "density" of the 
distribution is in this case unbounded. However, there is an old method 
[19] (see also [7]) that gives partial results. We should also mention that 
the Bernoulli problem has been solved in the continuum [5], but with 
methods that do not extend to lattice operators. 

In order to present the main ideas of our analysis, we introduce the 
function Na{E,uj) which counts the number of eigenvalues of less 
than E (recall that the Hamiltonian studied in this section is given by 
(10)). In the thermodynamic limit, an ergodic theorem guarantees the 
existence of the so-called integrated density of states 

N(E)= lim 

asz 3 \A\ 

which is independent of uj almost surely. Its derivative in E is the den- 
sity of states and is denoted by p{E). We can express Na{E,uj) by the 
resolvent Ga{E,uj) by noting that 

N A (E,u) =tr P A (E,u) = , P A (E, u)S x ) , 

x£A 
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where Pa(E, lu) denotes the spectral projection of /lyi(oi) onto the interval 
(—00, E), and 

1 r E 

Pa(E,u) =s-lim- / lm G A (E' +ie,u)dE' . 

e->0 TV J_ oc 

Using that 



2 r l 

* J-c 



dE 'l Zso 1 9 = ~ arctan(^-^) + 1 > @(E - x) , Ve > , 

, (x — E'Y + e z it e 



we see that the spectral theorem provides the bound 

N(E) = E (S xo , P(E, u)S X0 } < - f dE' Im E G(E' + ie; x ,x ), Ve > , 

*" J-oc 

(14) 

the first equality following from translation invariance. We can bound the 
probability that an eigenvalue of Via exists below E in terms of N(E), 

P[N A (E) > 1] = E1 {JV/1(E) > 1} < E^(S)1 {JVA(E) > 1} < ENa(E), 

and since Dirichlet boundary conditions raise eigenvalues we have that 

ENa(E) < \A\N(E) . 

The main idea is now to show that the integrated density of states in the 
band tails is so small that with very high probability, Ha does not have an 
eigenvalue below E. Then we can invoke the following well-known lemma 
to prove exponential decay of the off-diagonal elements. 
Lemma (Combes- Thomas). Whenever Ha has no spectrum below E\ then 
for E < E x , 

\Ga(E; x,y)\ < -exp(-constv / 5|x - y\), (15) 


where S := \E — E\\. 

Recalling equation (14), we see that, in order to get a bound on the 
integrated density of states, we have to bound the imaginary part of the 
averaged Green function. The strategy is simple. If v is small the first 
thought is to expand the resolvent (—A + v(uj) — E — ie) -1 in powers of 11 
around (—A — E) -1 . But since we consider the average E (—A — v(lu) — 
E)^ 1 , there may be an optimal energy Eo around which to expand. This 
energy is found as follows. We introduce a perturbation parameter A, 
writing 

h = -A + \v , 

with v — {vj}j eZl 3 a collection of independent Bernoulli random variables, 
Vj = ±1 with probability i, each; (for simplicity, we consider the symmet- 
ric case, x = |; but the general case is hardly more difficult). Formally, 
we have that 

G(E + ie) = - = - 

V ; -A + Xv-E-ie -A + E - is + Xv + (-E - E) 

= GoJ2((- Xv + E + E ") Gt >) n > 
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where we have introduced the unperturbed Green function 
G (-E + ie) = (-A + E - ie)" 1 . 

The first terms are 

G = Go - G [Xv -(E + E )]Go+ 

+ Go[\v -(E + E )]Gq[\v -(E + E )]Go + ... 

Using that Ev = and Ev 2 — 1, we obtain 

E G = Go + Gl[E + E ] + A 2 G (0, 0)G ( 2 , + G 3 [E + E } 2 + ... . 

Thus, we see that, in order for the A 2 -term to vanish, we must choose 

E + E = -X 2 Go(-E + ie; 0, 0) = 0(X 2 ) (16) 

to arrive at E G = Go + 0(\ 4 ). 

To make the above considerations mathematically respectable, we it- 
erate the second resolvent identity 

G(E + ie; x ,x) = Go(—E + ie; xo, x)+ 

J2 Go(-E + ie; x , y)[-\v + (E + E )]G(E + ie; y, x) , 

y 

with < E = -E - A 2 S(S), where A 2 S(S) is motivated by (16) and is 
defined in a self-consistent way by 

A 2 £(£0 = A 2 (-A -E- \ 2 T,(E) - ie) _1 (0, 0) . 

Since Eq has to be positive in order for Go(-Eo) to exist we have as an 
upper limit for E 

E* =A 2 A" 1 (0,0) <0. (17) 

Next, we iterate the resolvent identity M times with the intention of opti- 
mizing the truncation parameter M, later on, to minimize the remainder 
term. Setting W := -Xv - A 2 S and G £ := G(E + ie) we get 

M 

° £ = G o i WG o] m + G o [WGV\ M WG £ . 

m=0 

The dangling factor G e in the remainder is estimated trivially by 1/e. 
The key observation is that the imaginary part of the first M + 1 terms on 
the right-hand side can be shown to be proportional to e, whereas we will 
indicate below how to prove that the remainder term multiplying G £ is of 
order e~' 2A - ) . Thus we choose e 2 = e~* 2A ' and get the estimate 

ImEG(£' + i£;xo,a;o) < e~3< 2A )~ 1/2 . 
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In computing the terms of the form 
E[G (\vG r}(xo,x) 
=E[ ^2 Go(x ,xi)\v xl Go(xi,X2) ■ ■ ■ Xv Xm G (xm,x)] , 



it is easiest to use a graphical representation: Go(x,y) corresponds to a 
line joining x and y, while the interaction v corresponds to a vertex. Av- 
eraging over the randomness yields terms represented by graphs obtained 
by fusing an even number (since expressions involving an odd number of v 
vanish upon averaging) of vertices at a time until none remains unpaired. 
Because 

1 = Ev 2n « 7 (n) (Et; 2 )™ = 7 (n) 
7 (n) = = number of full pair contractions 

we can use Wick's theorem to get an upper bound by considering only 
fusions of pairs of vertices. It is well known how to bound such graphs: 
Consider the bubble graph: 



E J (VK> + |x|) 2 " 

Strings of n such bubbles are of size A n and add up to A/(l — A) (con- 
vergence only holds if Eo > A 4-2 * 5 ). 
It turns out that a graph of size M is of the order of 

M 

, up to logarithmic corrections. 



There are less than ~ 2 M (M/e) M ^2 < 2 M M! graphs generated by 

contracting pairs of 2M vertices, so the sum of all graphs is bounded by 



(4= 

\Ve7i 



M 

< M\(2X) MS . 



We now optimize our choice of the truncation parameter M: Choosing 
M = {2\y s we get the bound 

M\M~ M < Ce~ M = Ce- (2A) ^ . 

By similar calculations, collecting our previous estimates and setting 5 = 
1/2, we obtain the following upper bound for the integrated density of 
states, 



N(E) < Ce-^ 2X) 



-1/2 



We have thus shown that with probability larger than 



\A\Cc-^ 2X) 1/2 (18) 
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there is off-diagonal exponential decay of the resolvent Ga(E), for E G 
[—A, —A 2 ], see (17). Therefore, for any initial scale Zo there is a A such 
that (13) holds with high probability. Having established condition (13) of 
the multiscale analysis, we are looking for a Wegner estimate. Because of 
the very singular nature of the Bernoulli potential, the inductive scheme 
devised in [4] does not work. In [19], the following trick was introduced: 
For each scale l n , we define fj,(J) — (J,i n (J) to be the expected number of 
eigenvalues of h,A, n in an interval J. Next, we introduce the classical form 
of the Wegner estimate, and show that it implies (12). Since 

P[\\G Al (E)\\ > e !l/2 ] = P[dist(E, a(h Al )) < e^], 

it is clear that (12) holds with probability pi = 1 — l~ k if we can show 
that 

P[dist(E, <7(hA,)) <k]< C\Ai\k, 1/2 . 

The following easy estimate shows how we have to proceed: 

P[dist(E, a(h Al )) < k] < EN Al {E + k) - EN Al (E - k) = m(E - k, E + k) . 

We see that (12) is fulfilled if we exclude a set of "singular" energies, 
and the following lemma shows that this set of energies has very small 
measure. 

Lemma. Let n be a measure on an interval I. Let S be the set of energies 
E for which the measure is singular at scale e, that is at which 

n(E-e,E + e) > e 1/2 . 

If \S\ denotes the Lebesgue measure of S then 

\S\<2^I)e 1/2 . 

Proof. An easy application of Fubini's theorem shows that 

e 1/2 \S\ < J dE Jd»(x)l lE _ eiE+e] (x) 

< j d^ j dEl [E _^ E+e] {x) = A t(/)2e. 

□ 

Appealing to multiscale analysis, we see that, at each scale Z„ ~ Z§ , 

i l2 n 

we have to exclude energies of measure Cexp(— ^Iq ). Thus the total 
measure of energies we might have to excise is 

\E™\ = \\JET C \ < jrCeM-^l 2 ") < 2Cexp(-i/|). 

n n=0 

Because the integrated density of states is so small, we can choose lo 
to be exponentially large in A -1 / 2 (but not larger because of the factor 
Ai oc l A in (18)), and hence the set of energies we have to excise is of order 
exp(-Iexp(lA- 1 / 2 )). 
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A Cartan's Lemma 



Lemma. Let f(z) be a function analytic in the disc {z : \z\ < eR}, 
|/(0)| = 1, and let r\ be an arbitrary small positive number. Then we have 
that 

\({x G [~R,R] : \ f(x)\ < 5}) < 30^RS l ° sA *f w , 

where M/(r) = max| z | =r \f(z)\ and A(-) denotes Lebesgue measure. 

This lemma says that a function analytic in a disc that is bounded 
away from zero at one point is not too close to zero at "most" points. 
Since Ga{E) is a matrix, we need a generalization to matrix- valued func- 
tions, which is due to Bourgain [4] . 

Lemma. Let F(x) be a real analytic self-adjoint N x N -matrix function 
of x £ Q = [a, b] n , satisfying the following conditions (with m <JC TV, 
B 1 ,B 2 ,B Z >\) 

(1) T(x) has an analytic extension Y(z) to z G D" :=jz£C: \z— < 
f \a-b\} n with 

||r(z)||<Bi z G D" 

(2) There is a subset J of {1,2, ... , iV} such that \J\ < m and for all 
z G D n 

II (R{i,...,n}\j^(z)R{i,...,n}\j) II < B2 , 

where Rs denotes coordinate restriction to S 

(3) For some oj' G we have 

\\T{u')- l \\<Bz 

Then 

H({x G fl : ||r(a;) _1 || > K}) < C||s||oon|a - b\ n K~^^^^ , 

where \x denotes a measure with bounded density g with respect to 
Lebesgue measure, and C, c > are constants. 

Assumption (2) is clearly fulfilled in our case of T = h — E, the subset 
J being the lattice sites of the union of the bad cubes, and assumptions 
(1) and (3) are self-explanatory. 
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